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ABSTRACT 

We present follow-up observations with the Sunyaev-Zel'dovich Array (SZA) of optically-confirmed galaxy 
clusters found in the equatorial survey region of the Atacama Cosmology Telescope (ACT): ACT-CL J0022- 
0036, ACT-CL J2051+0057, and ACT-CL J2337+0016. ACT-CL J0022-0036 is a newly-discovered, massive 
O 10 15 M ), high-redshift (z = 0.81) cluster revealed by ACT through the Sunyaev-Zel'dovich effect (SZE). 
Deep, targeted observations with the SZA allow us to probe a broader range of cluster spatial scales, better 
disentangle cluster decrements from radio point source emission, and derive more robust integrated SZE flux 
and mass estimates than we can with ACT data alone. For the two clusters we detect with the SZA we compute 
integrated SZE signal and derive masses from the SZA data only. ACT-CL J2337+0016, also known as Abell 
2631, has archival Chandra data that allow an additional X-ray-based mass estimate. Optical richness is also 
used to estimate cluster masses and shows good agreement with the SZE and X-ray-based estimates. Based on 
the point sources detected by the SZA in these three cluster fields and an extrapolation to ACT's frequency, we 
estimate that point sources could be contaminating the SZE decrement at the < 20% level for some fraction of 
clusters. 

Subject headings: X-rays: galaxies: clusters: general - cosmology: observations - cosmic microwave back- 
ground 
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1. INTRODUCTION 

The Sunyaev-Zel'dovich effect (SZE) is a small (typically 
< 1 mK) distortion of the cosmic microwave background 
(CMB) spectrum caused by the inverse Compton scattering 
of CMB photons off energetic electrons of the hot intracluster 
medium of galaxy clusters (Zeldovich & Sunyaev 1969; Sun 



yaev & Zerdovich||1970| |1972| for reviews see, |Sunyaev & 



Zeldovich 1980, Rephaeli [T995] |Birkinshaw| 1 999jjCarlstrom 
et al.|2002|>. The redshift independence of the SZE makes it 



a potentially powerful tool with which to search for galaxy 
clusters, especially in the distant universe. Abundances of 
clusters probe the growth of structure and have placed use- 
ful constraints on the fluctuation amplit ude, erg, and the mat- 
ter density, (e.g.,|Henry & Arnaud|1991||Viana & Liddle 
T9%l[B ahcall et alJ1997jlEke et al.|1998Por gani et al .|2001[ 
Reiprich & Bohringer 2002; Schuecker et al. 2003, Henry 
2004 



Mantz et al.||2008| |Vikhlinin et al.||2009 



Rozo et al.|2blOI |Vanderlinde et al.|20T0 
Tot cli 



Mantz et al. 



Sehgal et al. 



2010 

201 1 l The evolution of cluster abundance with redshift is one 
of the few probes of the growth of structure and has the po- 
tential to tightly constrain cosmological parameters and pro- 
vide insight into the equation of state of the dark energy (e.g., 



Bartlett & Silk 1994 ;|Holder et al.|2000l|Haiman et al.|200T| 



Majumdar & Mohr 2004|. The biggest challenge to realiz- 
ing the cosmological potential of cluster surveys is relating 
cluster mass to an observable such as integrated SZE signal. 
Well-determined masses for even a subsample of clusters can 



significantly improve parameter constraints (e.g., Majumdar 
|& Mohr|2004l 

The last few years have seen significant advances in surveys 
using the SZE. The Atacama Cosmology Telescope (ACT; 
Fowler et aL]p0a7] l, the South Pole Telesc ope (SPT; |CarL| 
strom et al.||2011| l, and the Planck satellite (Planck Collab- 



oration ||2011a[ l are, for the first time, produc ing catalogs 
of galaxy clusters discovered through the SZE (Vanderlinde 



et al |2010H Marriage et al. 2010 1 Planck Collaboration 201 lb ; 
Williamso n et al.||2011| l. The fast mapping speeds and im- 
proved sensitivities of these instruments have enabled them to 
survey large regions of the sky with sufficient depth to detect 
massive clusters, but their limited angular resolutions (> 1') 
do not allow detailed studies of cluster astrophysics with their 
data alone. 

We have performed initial follow-up observations with 
the Sunyaev-Zel'dovich Array (SZA) of three optically- 
confirmed ACT clusters in ACT's equatorial strip. Deep, tar- 
geted SZA observations of ACT clusters provide higher sen- 
sitivity over a broader range of cluster spatial scales than that 
of ACT. In addition, the spatial filtering of the interferometer 
provides a method of cleanly disentangling radio point source 
emission from the cluster SZE signal. ACT observational 
details and target selection are presented in Section [2] SZA 
and Chandra observations and data reduction are described in 



Sections 3.1 and 3.3 respectively. The analysis method and 
results are reported in Section|4] the implications in Section|5j 
and conclusions in Section|6] 

Throughout this paper, all uncertainties are reported at 68% 
confidence and we adopt a flat, A-dominated cosmology with 
Q M = 0.3, Q A = 0.7, and H = 70 km s _1 Mpc" 1 consistent 
with re cent Wilkinson Microwave A nisotropy Probe (WMAP) 
results ( |Komatsu et al.|201 l||2009] l. 



The Atacama Cosmology Telescope (ACT) is a 6-m diame- 
ter telescope located at an elevation of 5200 m in the Atacama 
desert of Chile. Three 1024-element transition-edge-sensing 
bolometer arrays operate at 148, 218, and 277 GHz and sur- 
vey large regions of the sky mainly in two regions, a south- 
ern strip centered at declination -52.5° and in an equatorial 
strip that encompasses th e Sloan Digital Sky S urvey (SDSS) 
Stripe 82 (hereafter S82, |Abazajian et al.|2009) . For the ACT 
instrument, observat i on, r eduction, and calib ration details se e 
IFowler et a l. ( 2007 ) , |Swetz et al.| ( |20TT| >, |Das et al.| ( |20TTT ), 
and |Hajian et al.| ( |2010[ ). Initial resul ts from the A CT south- 
ern strip in clude CMB powe r spectra (|Fowler et al.|2010{ |Das| 
|et al.|20lT|), compact source ( Marriage e t al.|201 1 J > and cluster 
catalogs (Marri age et al.|2010| l, cluster follow-up (Menanteau 
|et al.|2010| l, and anal ysis oFthe cosmolo gical implications of 
CMB power spectra (Dunkley et al.|20T0] > and cluster yields 

( |Sehgaletall20TT] i. 

A ma tched filt er method (|Haehnelt & Tegmark|1996 Her- 



|ranz et aL|2002| |Melin et al.|2006| l is used for cluster detec- 
tion to produce a cluster catalog (for ACT-specific details, see 



2. ACT TARGET SELECTION 



Marriage et al. 2010) . We focus on the S82 region as op- 
tical data are available from which we measure cluster red- 
shifts and estimate cluster masses via richness and luminos- 
ity. We chose three preliminary high signal-to-noise (S/N* 5) 
cluster detections from the S82 region for follow-up as a pi- 
lot study for SZA observations; these were ACT-CL J0022- 
0036, ACT-CL J205 1+0057, and ACT-CL J2337+0016. We 
will refer to these clusters as ACTJ0022, ACTJ2051, and 
A2631, where ACT-CL J2337+0016 is a previously known 
Abell cluster. The upper row of panels of Figure [TJ shows the 
SZE images of these three fields expressed as Compton-y (see 
Eq.|2]i made from recent ACT maps. ACTJ0022 is a new high- 
redshift cluster found by ACT through the SZE. ACTJ0022 
was spectroscopically confirmed at z — 0.81 from long slit 
observations of the brightest cluster galaxy (BCG) conducted 
at the Apache Point and Gemini South Observatories (more 
details on these observations will be described in Menanteau 
et al. 201 1, in prep., where we present the X-ray and optical 
properties of the S82 ACT clusters). The S82 data provide the 
spectroscopic redshift, 0.33, fo r ACTJ205 1 and A263 1 has a 
spectroscopic redshift of 0.27 (Str uble & Rood|[T 999). New 
maps produced with the full ACT equatorial data set reveal 
ACTJ2051 to be lower significance (S/N 4) than inferred 
from the preliminary maps. Table [T] lists the clusters, alter- 
nate names, and redshifts. 



3. DATA 

3.1. SZA Observations 

The Sunyaev-Zel'dovich Array (SZA) is an 8-element in- 
terferometer with 3.5-m diameter dishes operating at 30 and 
90 GHz. It has an 8 GHz bandwidth, and is located at an alti- 
tude of 2200 m in the Inyo Mountains of California. Our ob- 
servations use the 30 GHz system, which has a 10.'5 FWHM 
primary beam (field of view). Typical system temperatures 
are T sys 40-50 K, including atmospheric contributions. Six 
of the antennas are placed in a compact configuration to max- 
imize cluster sensitivity, while two are deployed at longer 
baselines to provide higher resolution data. The spatial filter- 
ing of the interferometer allows one to disentangle the small- 
scale, positive point source emission from the large-scale, 
negative SZE signal at these freque ncies (for an example in 
a similar context see Reese et al. 2002 ). The SZA has become 
part of the Combined Array for Research in Millimeter-wave 
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TABLE 1 
SZE Image Statistics 
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• Pre- and post-CARMA transition of the integration of SZA with CARMA. 
b Effective on-source integration time accounting for excised data. 



Astronomy (CARMA) and our data were taken both before 
and after the CARMA transition. 

The SZA observed the three ACT clusters between March 
and May and between August and September 2010. A bright 
quasar was observed every 30 minutes, for about 2 minutes 
integration time, to monitor the system gain and phase during 
each observation. Since the fluxes of quasars can be variable 
on timescales of days or months, these sources only serve as 
sec onda ry calibrators. We discuss absolute calibration in Sec- 
tion [321 

Pre-CARMA data are reduced using a suite of MATLAB 
routines dev eloped by the SZA collaboration (e.g., Muchovej 
et al.||2007|). Post -CARMA data are reduced using Miriad 
Data are excised when one telescope 



(Sault et a 



1995! 



shadows another, when cluster field data are not straddled by 
two phase calibrator observations, when there are anomalous 
changes in instrumental response between calibrator obser- 
vations, when the system temperature changes dramatically 
between integrations, or when there is spurious correlation. 

The Fourier plane is also known as the u-v plane, where 
(m, v) are the Fourier conjugate variables to right ascension 
and declination, respectively. We define the radius in the u-v 
plane as r uv - Vm 2 + v 2 , which corresponds to the baseline 
length in units of the observation wavelength, A. The con- 
figurations used for these observations produce a break in u-v 
coverage at 2 kA, with similar sensitivities for the low (r uv < 2 
kA) and high (r uv > 2 kA) resolution data (15 short and 13 
long baselines). Table[T|summarizes the on-source integration 
times, and (for the high and low-resolution data separately) 
sensitivities and the FWHMs of the synthesized beams (effec- 
tive resolution). The rms map noise, <x CMB , in CMB tempera- 
ture units for the given synthesized beams, is also given for the 
low-resolution maps. The total integration times reported cor- 
respond to the equivalent time that the full 8 -element, 8 GHz 
bandwidth SZA spent on each cluster (accounting for excised 
data). 

Data visualization is done with Difmap (Shepherd 1997|l . 
Identification of point sources in the field is performed with 
the high-resolution (r,„. > 2 kA) data. Point sources are mod- 
eled to remove them from the data, and the low-resolution 
data (r UY < 2 kA) are then deconvolved to produce images of 
the cluster fields with sources subtracted. We present decon- 
volved images of the point source-removed SZA observations 
as a measure of the data quality. The model fitting is per- 
formed in the Fourier plane directly. The cluster and any point 
sources in the field are modeled simultaneously using the data 
in its entirety. Mo deli ng of the cluster and point sources is dis- 
cussed in Section P+TTl and results are discussed in Section l4~2l 

The lower panels in Figure [T] show the deconvolved SZA 
images of the three cluster fields after removal of the point 



TABLE 2 

Pre- and post-CARMA Point Source Flux Densities 











Field 


(mJy) 


(mJy) 


Ratio 


ACTJ0022 


0.95 ± 0.20 


0.80 ±0.14 


1.19 ±0.30 




1.07 ±0.18 


0.65 ±0.13 


1.65 ±0.36 


ACTJ2051 


8.88 ±0.33 


5.03 ±0.17 


1.77 ±0.32 




3.66 ± 1.36 


2.49 ± 0.61 


1.47 ± 1.05 


Weighted Avera 






1.51 ±0.07 



sources. Contours are multiples of twice the rms noise in 
each map, and the color scale is in units of Compton-}'. 
The FWHM of the synthesized beam (effective resolution) is 
shown in the lower left of each panel. The images all cover 
the same angular scale (1 1' on a side) and have the same color 
scale to facilitate comparison. The ACTJ0022 and A2631 
fields reveal high-significance SZA detections. The point 
sources in the ACTJ205 1 field are used both to strengthen our 
calibration of the SZA data and to better assess SZE decre- 
ment contamination due to such sources. 

3.2. SZA Calibration 

Pre-CARMA SZA operation s routinely m onitored Mars for 
absolute calibration (e.g., [Muchovej et al. 2007 1 



In the pe- 
riod right after the SZA-CARMA merger, standard calibration 
protocols had not yet been implemented. Therefore, though 
the system gains are stable, the absolute calibration is off by 
an arbitrary multiplicative factor. We calibrate the transitional 
post-CARMA data using radio sources observed both before 
and after the integration to derive an average calibration fac- 
tor that is applied to the data. This is non-ideal because the 
sources can vary; however, we note that on average we do 
not expect the source flux densities to exhibit any particular 
trend upwards or downwards over time, so the average ratio 
of several sources pre- and post-CARMA may be expected to 
follow any changes in calibration. The cluster signal will re- 
main constant, so we can validate our calibration by fitting the 
same model to the pre- and post-CARMA data and comparing 
the central Compton-y values. 

The ACTJ0022 and ACTJ2051 fields have both pre- and 
post-CARMA data. For the calibration analysis, point source 
models are fit to the data. A cluster model is simultaneously 
fit to the ACTJ0022 data. A Markov Chain Monte Carlo 
(MCMC) analysis is performed to derive flux densities and 
unc ertain ties (model and fitting details are described in Sec- 
tion |4Tj. Ratios of pre- to post-CARMA flux densities are 
computed. The flux densities and ratios are summarized in Ta- 
ble[2] The average ratio is 1.5 ±0.1, computed with an inverse 
variance weighting. There is a wide source-to-source varia- 
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Fig. 1. — Upper. ACT Compton-y images of the three optically confirmed clusters made from recent ACT maps. We fit (excluding the cluster region in the fit) 
and remove a quadratic polynomial from the data, smooth with a Gaussian, and resample to smaller pixels. The contours are multiples of twice the rms of each 
map (black). Also shown are contours from the match filtered map (gray) as multiples of twice the rms. Lower: Deconvolved SZA images of the three cluster 
fields in units of Compton-y. The SZA images are made with short baseline data < 2 k/i after removal of point sources. Contours are multiples of twice the rms 
of each map (for details of the SZA observations presented here, see Table[l). The color scales are in units of 10 5 Compton-y, and solid (dashed) lines represent 
positive (negative) signal. All figures cover the same angular scale (11' on a side) and are on the same color scale to facilitate comparison. The FWHM of the 
effective resolution (synthesized beam for the SZA) is shown in the lower left corner of each panel. For visualization purposes the SZA images incorporate only 
the low-resolution SZA data and the effective beams reflect this choice. The SZA data analysis (Sectionpl uses the data in its entirety. ACTJ0022 and A2631 
both yield high S/N detections. ACTJ2051 is not detected by the SZA, which is consistent with its lower SZE signal in the current ACT maps (shown here) and 
with that expected from optical data (see Section[5j. In general, the qualitative agreement between the ACT and SZA data is good. 



tion in the ratios, which range between 1 .2 and 2.0, though the 
uncertainties on the ratios are relatively large and the largest 
discrepancy is < 2cr. Our average ratio is consistent with the 
rescaling 1 .4) found by members of the CARMA collabo- 
ration (T. Plagge, private communication). 

The cluster SZE signal will not change, but comparisons are 
complicated by its extended structure. We also fit the same f3- 
model (fixed core radius and power law index) to the pre- and 
post-CARMA data for ACTJ0022. The ratio of the central 
Compton-y parameters also yields the same factor of 1.5. That 
the cluster normalizations yield the same calibration scaling 
as the point sources lends some support to our simple ap- 
proach. As a result, the post-CARMA data are simply scaled 
by 1.5 and used for the analysis. Note that the rms values in 
Tablefllcorrespond to the scaled post-CARMA data. The pre- 
CAR MA data absolute flux c alibration is known to better than 
10% (Muchovej et al.||2.007) >. Given the large scatter in the 
flux density ratios, we adopt a conservative 20% calibration 
uncertainty to account for potential calibration systematics. 



3.3. Chandra Observations of A2631 

ACTJ2337 is the Abell cluster A2631, which has publicly 
available Chandra X-ray observations. The data consist of 
two ACIS-I observations, obsIDs 3248 and 11728, of dura- 
tions 9 and 17 ks respectively. The data are reduced start- 
ing with the level 1 events file using CIAO 4.3 and calibra- 
tion database version 4.4.1. Standard corrections are applied, 
along with light curve filt ering and other stan dard processing 
(for reduction details see Reese et al. 2010). The Chandra 
images of A2631 are made by binning the 0.7-7.0 keV data 
to 1'.'968 pixels and exposure maps are computed at 1 keV. 
The images and exposure maps from the two observations are 
combined, and a wavelet based source detector is used on the 
combined image and exposure map to identify and generate a 
list of potential point sources. The list is used as the basis of 
our point source mask. 

Figure [2] shows the smoothed Chandra data (color) along 
with the deconvolved SZA data (contours). The color scale 
shows Chandra counts and the contours are multiples of twice 
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Fig. 2. — Deconvolved SZA image of A263 1 (contours) on the smoothed 
Chandra counts image (color). Contours are multiples of twice the rms of 
the map (see Table u\ for details) and the color scale is in units of Chandra 
counts. The FWHMof the synthesized beam (effective resolution) is shown 
in the lower left corner. The E-W elongation is seen in both the SZA and 
Chandra data, and the peaks in the SZE and X-ray data align well. 

the rms of the SZA map. The FWHM of the synthesized beam 
(effective resolution) is shown in the lower left corner. The 
same E-W elongation seen in the SZA data is seen in the 
Chandra data as well, and the peaks in the SZE and X-ray 
images align. 

4. ANALYSIS 
4.1. Method 

We fit cluster models to the data and derive parameters from 
the best-fit models. ACTJ0022 has no X-ray data and we fit 
models to the SZE data only. A263 1 has both SZA and Chan- 
dra data and we perform three complementary analyses using 
1) SZA data only; 2) Chandra data only; and 3) both the SZA 
and Chandra data jointly. Derived quantities such as mass and 
integrated SZE signal are computed within Z?5oo, the radius 
within which the mean interior density is 500 times the critical 
density at the redshift of the cluster, p c (z) = 3i/ 2 (z)/(8^G), 
where H(z) is the Hubble parameter, and G is Newton's grav- 
itational constant. 

4. 1 . 1 . Chandra Spectroscopy for A2631 

Both a single temperature and a temperature profile are 
measured for A2631. We extract spectra and response files 
for both data sets separately. Single temperature spectra are 
extracted within (0.15-l)7? 5 oo for the isothermal temperature 
analyses. The region used for spectroscopic extraction is 
found recursively, by picking a trial radius, extracting a spec- 
trum to determine the electron temperature, T e , within that 
radius, and computing M500 = M(R^o) from hydrostatic equi- 
librium (Eq. [3]) and the resulting /?5oo- This process is re- 
peated until the values of the input and output R^oo agree to 
1%. Spectra are also extracted from within the full Rsqq deter- 
mined as outlined above for comparison. 

For the temperature profile, we extract spectra in rad ial an- 
nuli with r out /Vj n set to a constant ratio. This follows [Vikhlinin| 



TABLE 3 
A263 1 Chandra Spectral 
Properties 
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Note. — Fits performed with 
redshift z = 0.27 (Struble &| 
|Rood|1999) and Galactic H I col- 
umn density N» = 3.55 x 10 20 
cm -2 dKalberla et al. 12005 1. 



et al. ( 2006| hereafter V06) who found this choice produces 
roughly equal counts per annulus for cluster observations. 
However, we choose r out /r; n = 2 (instead of 1.5 used by V06) 
in order to increase the counts in each radial bin, since our 
data have many fewer counts than the data considered in V06. 

We construct quiescent ACIS background datasets by re- 
projecting the ACIS blank-field observations to match each 
dataset. To account for variations in the particle background, 
the blank-field observations are scaled by the ratio of fluxes 
in the 9.5-12 keV band, where the Chandra effective area is 
essent ially zero and the flu x is due entirely to the background 
(e.g., Vikhlinin et al.|2005 ). Background spectra are extracted 
from these quiescent background data sets using the same re- 
gions as that for the cluster data and used when fitting the 
cluster spectral properties. 

Spectra for both Chandra observations using their respec- 
tive responses are fit simultaneously to the same plasma 
model with the normalizations allowed to vary independently 
betwe en data sets. Xspec (Arnaud 1996] [Dorman & Arnaud 
200 l|l is used to model the intracluster medium with an APEC 

that includes Bremsstrahlung 



spectrum (|Smith et al. 2001 



We a dopt the Galactic col - 



3.55 x 10 20 cm' 2 dKalberla et al.|2005b 



and line emission components 
umn density of 7V H 
so lar abundances of fA splund et al. ( 200 9|), an d cross sections 
of|Balucinska-Church & McCammon ( 1992) with an updated 

|1998|l. The analysis uses data in 
; 'cstat" statistic, which is 



He cross section ( Yan et al.[ 
the 0.7-7.0 keV energy range. The 

similar to the |Cash| (|T979 ) statistic, is used when modeling the 
spectral data to properly account for low counts. 

Table|3]summarizes the Chandra spectral results, listing the 
inner andouter extraction radii, derived temperatures, T e , and 
metallicity relative to solar, Z. In all cases, the model provides 
a good fit to the data, without any obvious structure or pattern 
in the residuals. As a check, the T e profile spectra are fit fixing 
the metallicity to the global value (Z = 0.4). The resulting 
TVs agree with those in Table[3]to within lcr for all regions. 

4.1.2. Cluster Models 

For the SZA data, the clusters are fit with both the tradi- 
tional y6-model (Cav aliere & Fusco-Femiano|1976||1978[ ) and 
the universal pressure profile of Arnaud et al.|(|2010 hereafter 
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A10). For the SZE-only jfi-model analy sis, we fix yg = 0.8 6 
from a best fit to an average SZE profile (Plagge et al.|2010 



For the Chandra data, the cluster is fit with the traditional 
isothermal yS-model and the modified /3-model (e.g., V06) for 
density along with the corresponding V06 T e (r) profile, here- 
after referred to as the V06 model. Many of the V06 model 
parameters must be fixed, given the large number of parame- 
ters compared with the quality of the X-ray data. A263 1 does 
not show signs of a cool core so we excise those parts of the 
V06 model. 

Joint fits to the SZA and Chandra data are performed by 
assuming the A10 pressure profile and a simplified, core-cut 
form of the V06 density profile (hereafter sV06) for the cluster 
models. These fits are performed both including and exclud- 
ing the X-ray spectroscopic data, allowing the SZE and X-ray 
surface brightness data to place constraints on the c luster tem- 
perature (for details see Mroczkow ski et al.|2009[ ). The tem- 
perature is derived assuming the ideal gas law (P e = n e k B T e ). 
The X-ray surface brightness data determine the sV06 density 
fit, and the SZE surface brightness data drive the A10 pressure 
profile fit. A single temperature i s derived from t he in ferred 
T e profile using the prescription of Mazzotta et al. (2004} over 
the (0.15-l)/?5oo region. For the AlO+sV06 fits that include 
spectroscopic information, the likelihood from the inferred T e 
described above is included in the MCMC using the output 
from the Chandra single-temperature spectroscopy. A sep- 
arate MCMC run is done without including the T e spectral 
component in order to assess how X-ray surface brightness 
and SZE data can be used to find M(r) without relying on 
X-ray spectroscopy. 

4.1.3. Radio Point Source Model 

The radio point source model accounts for the primary 
beam attenuation and includes a spectral index that models 
the frequency dependence. The spectral dependence of the 
point source model is given by 



F v — F v, 



(i) 



where F v is the intrinsic point source flux density at frequency 
v, F Vo is the intrinsic flux density at fiducial frequency vo, and 
a is the spectral index. We adopt the average frequency for the 
fiducial frequency, vo = 30.938 GHz, and report flux densities 
at this frequency. 

Point sources are first modeled with both F Va and a allowed 
to vary. The 8 GHz bandwidth of the SZA potentially pro- 
vides leverage on both the spectral index and flux density of 
each point source. The flux density is, in general, well con- 
strained. However, only weak constraints are obtained on a 
for all but the brightest sources. We first run fits including a in 
the Markov chain and then determine a using our 30 GHz flux 
density combined with 1 .4 GHz flux densities fr om the NRAO 
VLA Sky Survey (NVSS; [Condon et al p998) . As a check, 
we compare the flux densities obtained from the NVSS-based 
spectral indices to those obtained just from the SZA data. The 
flux densities are consistent within the 68% confidence re- 
gions. 

4.1.4. Markov Chain Monte Carlo Analysis 

We perform a Markov chain Monte Carlo (M CMC) anal- 
ysis of the SZE and X-ray da ta (for details see Reese et al. 
2000 Bonamente et al. 2004). The philosophy behind the 
analysis is to keep the data in a reduced but native state, and 



to run the models through the observing strategy to compare 
directly to the data. Interferometers measure the Fourier trans- 
form of the sky brightness modulated by the primary beam. 
Therefore the SZA data most naturally provide constraints in 
the u-v plane, where the noise properties of the data and the 
spatial filtering of the interferometer are well understood. We 
perform model fits directly in the u-v plane, computing the 
Gaussian SZE likelihood, Xsze- For X-ray data, the Poisson 
likelihood, _£x, is computed for each pixel, ignoring those that 
fall within the point source mask. The SZE and X-ray data are 
independent, and their likelihoods can simply be multiplied 
to obtain the combined likelihood. Best-fit parameters and 
confidence intervals are determined from the 50%, 16%, and 
84% levels of the cumulative distribution function, which de- 
fine the the median and 68% confidence region. The resultant 
probability distributions are visually inspected, and conver- 
gence and mixing of the Markov c hains are checked with the 
Geweke Z-statistic (Gew eke)l992| l. 

The cluster and any detected point sources in the field are 
modeled for the SZA data. The SZE signal varies as 



Ar szE = f(v)T c] 



dl- 



k K T„ 



' T n e o- T = f(v)T CMB y, (2) 



where /(v) is the frequency dependence of the SZE at fre- 
quency v, n e and T e are the electron number density and tem- 
perature, k B is the Boltzmann constant, m e is the mass of the 
electron, <x x is the Thomson cross section, integration is along 
the line of sight (, and y is the Compton-y parameter. SZA 
data have sixteen 500 MHz bands covering 8 GHz of band- 
width. The frequency dependence of the SZE is taken into 
account when modeling. Relativistic cor rections to f(v) (e.g., 
Itoh et al.|1998)|Challinor & Lasenby[T998] l depend on T e and 
are not included in this analysis for consistency, since only 
one cluster has a measured T e . The effects of this small (< 3% 
at 30 GHz) correction are discussed in Section[5] 

For A2631, Chandra data are modeled with a cluster model 
and a constant X-ray background. Regions containing point 
sources are excluded from the fit. 

4.1.5. Derived Cluster Properties 

Derived quantities such as the integrated Compton-y, isoo, 
and mass, M500, are computed for each step in the Markov 
chain for each of the above types of fits we perform. We 
use these output chains to determine the median values and 
68% confidence regions for each parameter of interest. This 
method cleanly propagates uncertainties from the parameters 
included in the chain. For example, uncertainties from mod- 
eling detected radio sources are propagated into the integrated 
Compton-y results. 

We compute the integrated Compton-y, T sph (r) = 
cr T /(m e c 2 ) J P e dV, within a spherical volume enclosed by ra- 
dius r, where dV = 4nr 2 dr and P e = n e k B T e is the elec- 
tron pressure. Since this latter quantity tracks thermal en- 
ergy (E = 3/2 JP eixi dV), and thermal pressure is the domi- 
nant source of su pport against gravitational collapse (see, e.g. 
Nagai et aL|2007| |Lau et al.||2009| who report that only 10- 
20% of the pressure comes from non-thermal support), we 
can expect Y spb {r) to track gravitational energy within radius r. 
Assuming the virial relation, a constant gas fraction of 0.13, 
an average metallicity, and a total mass profile that can be de- 
scribed by a Navaro, Frenk , and White (NFW) halo model 
( |Navarro et al.|[T996l |1997| l, we estimate mass from fits to 
Y sph (r) derived from our A10 and j6-model fits to the SZE 
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TABLE 4 
Radio Point Sources 



Field 


RA. (J2000) 
(h m s) 


Decl. (J2000) 
(d m s) 


R" 
(') 


(mJy) 


(mJy) 


a 


pest 
1 148 

(mJy) 


ACTJ0022 


00 22 13.006 
00 22 05.987 


-00 36 33.35 
-00 35 54.95 


0.3 
1.6 


1 33+° 15 
98+° 13 

-0.13 


18.8*? 

33 2 +L1 
—l.l 


-0.9 ±0.1 
-1.1 ±0.1 


0.3 
0.2 


ACTJ2051 


20 51 17.758 
20 51 40.720 


+00 53 20.16 
+00 52 02.74 


2.3 
7.3 


7 36+°- 14 
4 01+ 044 


43.8+I; 4 
! 6 9+0.7 


-0.6 ±0.1 
-0.5 ±0.1 


3.0 
1.9 


A2631 


23 37 40.209 
23 37 38.106 


+00 16 42.18 
+00 10 01.93 


0.7 
6.2 


T. qq+0.16 
3 20+ 034 

J - -0.33 


107.0+ 3 jj 

i j 4 +0.6 
^ -0.6 


-1.1 ±0.1 
-0.4 ±0.1 


0.8 
1.6 



' Projected distance from the cluster center. 



TABLE 5 
SZE and X-ray Derived Properties 



Field 


^500 


M 5m 


Ysoo 


(model) 


(Mpc) 


(10 14 M G ) 


(If)- 4 Mpc 2 ) 


ACTJ0022 








/5-SZE 


1 07+ 004 

1,u '-0.04 


8.8!|; 


1.04*3° 


A10-SZE 


1 00+ 005 

'•""-0.05 


7.3!i;0 


0.84+0;!$ 


A2631 








/3-SZE 


1 42+ - 15 

1/tZ -0.13 


11.5+ 4 ' 


i aa+0-85 


A10-SZE 


i co+0.36 
1.58_ 025 




i co+1.41 
1.58_ 065 


/J-X-ray 


i -3 1+0.03 
J -0.03 


a 5+O.8 




V06-X-ray 


1 54+O.IO 


13.8+1' 




A10+sV06 a 


i oo+0.20 
-0.12 


Q 4+4.8 


1 03+ ' 3 ' 


A10+sV06 b 


1 31+ - 35 

L,D -0.38 


9 1+ 9 - 4 


1 01+ 065 



" Includes X-ray spectroscopic information. 
b Does not include X-ray spectroscopic informa- 
tion. 



data (as was done in Mroczkowski 2011). These SZE-only 
mass estimates can be performed regardless of chosen model 
fit, and, like estimates from the X-ray assuming hydrostatic 
equilibrium, rely on spherical symmetry and do not take into 
account sources of non-thermal pressure support. We define 
F500 = T spll (,Rgoo) as the integrated Compton-3' within R^qq. 

Mass estimates from th e Chandra data are based on hydro- 
static equilibrium (e.g., S arazin|1988 1, 

r dP Jr) 



M(r) = -- 



k B T e (r)r ( dln(n e ) dln(T e ) 



G/irrip \ d\n(r) d\n(r) 



Gp gas (r) dr 



where /i is the mean molecular weight, m p is the mass of the 
proton, and p gas and P gas are the total gas density and pres- 
sure, respectively. The mass as a function of radius is used to 
compute R 500 using M 500 = M(R 500 ) = (4jr/3)R i 5Q0 500p c (z), 
where p c (z) is the critica l density at red shift z. Our analysis 
is similar to that of Vik hlinin et al.| ( |2006l >. 



4.2. Results 

The MCMC results for all detected point sources are pre- 
sented in Table |4] In addition to the SZA 3 1 GHz flux densi- 



ties, the corresponding NVSS ( [Condon et al.|1998| > 1.4 GHz 
flux densities are listed. These two flux densities are used 
to compute the spectral index a, which is then used to esti- 
mate the point source flux density in ACT's 148 GHz band, 
Fm (both also listed in TableBl. The projected distance from 
the cluster, R, is also listed and gives an idea of the potential 
impact of each source on cluster detection and potential con- 



tamination of the SZE flux. Table [5] summarizes the cluster 
modeling results for the SZA, Chandra, and joint analyses. 
We report R500, and T500 from our MCMC analysis. 

Figure [3] shows the SZE u-v radial profiles of the SZA 
data (points) with best-fit /3 (blue) and A10 (red) models 
for ACTJ0022 (left) and A2631 (right). The visibilities (u- 
v plane data) are converted to a frequency independent u-v 
plane Compton-y, Y uv , and scaled by the angular diameter 
distance squared, D\, cr eating a SZE lumin osity-like quan- 
tity (for details see, e.g., Mroczko wski et al.|2009| l. The real 
parts of Y uv are plotted. The imaginary components are con- 
sistent with zero. Residuals are shown in the lower sections 
of both panels. The ACTJ0022 radial profile has x 2 = 35.1 
and 35.8 with 28 degrees of freedom (dof) for the /3 and A10 
models, respectively. The corresponding probabilities of ob- 
taining these x 2, s or larger by chance given the degrees of 
freedom are, p(> x 2 \ dof) = 0.17 and 0.15. The A2631 radial 
profiles have X 2 = 38.1 and 36.1 with 28 dof for the/? and A10 
models, respectively, with corresponding p(> x 2 \ dof) = 0.10 
and 0.14. 

Figure|4]shows the Chandra surface brightness profile (up- 
per) and temperature profile (lower) along with the best fit 
models. The data (points) and best fit /3 (blue) and V06 (red) 
models are shown. The best-fit background level is shown by 
the dotted line. The lower portion of each panel shows the 
residuals in units of the standard deviation. Despite this clus- 
ter's E-W elongation (Fig. J2l, a spherical model, on average, 
provides a good fit with^ = 190.3 and 155.9 for the /J and 
V06 radial profiles that have 194 and 191 degrees of freedom 
(dof), giving p{> X 2 \ dof) = 0.56 and 0.97, respectively. The 
V06 temperature profile results in x 2 = 2.6 for 2 degrees of 



(3) freedom with p(> x 2 \ dof) = 0.27. 



4.3. Optically Informed Cluster Properties 

SDSS S82 data provide redshift information and enable 
mass estimates based on relations between cluster ma ss and 
optical properties such as richness and luminosity (e .g., |John-| 
|ston et al. 120071 |Reyes et al.|2008||Rozo et al.|2009) . We have 
computed cluster masses from the A^oop-^soo scaling relation 
of Rozo et al. ( 2009 ), which is based on weak-lensing and X- 



ray measurements of clusters from the SDSS maxBCG cata- 
log ( Koester et al.|2007 1. Masses computed from the A^2oop- 
M500 scaling relation of Reyes et al. (2008) yield similar re- 
sults after correcting for the different radii used between the 
studies. M500 is the mass within R500, as defined earlier, the 
radius within which the mean interior density is 500 times that 
of the critical density at the redshift of the cluster. A^200p is the 
number of red sequence galaxies in a cluster measured within 
RiOQp, the radius within which the mean interior density is 200 
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Fig. 3. — SZE profiles for ACTJ0022 (left) and A2631 (right) as a function of u-v radius, r m = Vn 2 + v 2 . Best-fit A10 (red) and /? (blue) model fits are also 
shown. Residuals are shown in the lower sections in units of the standard deviation. Point sources have been subtracted from the visibilities and the phase centers 
shifted t o the cluster centers before computing the profiles. The u-v plane Compton-y, Y uv , is scaled by D\, creating a SZE luminosity-like quantity (for details 
see, e.g., Mroczkowski et al. 2009 1. The real parts of F,„, are shown. The imaginary parts are consistent with zero. The dotted line shows the shadowing limit of 
the 3.5 m dishes. 



Radius (kpc) 
100 
T" 




Radius (arcmin) 



Fig. 4. — Upper: A263 1 Chandra surface brightness profile (points) with 
best fit fi (blue) and V06 (red) models. The lower section shows residuals in 
units of the standard deviation. The models are similar and the simple spheri- 
cal /3-model provides a good fit to the Chandra data. Lower: A263 1 Chandra 
temperature profile (points) with best fit V06 temperature (red) model. The 
lower section shows residuals in units of the standard deviation. 



times the mean matter density at the redshift of the cluster, 
p(z) = Q.m{z) p c {z), and is denoted by the subscript p. T^oop" 
at z — corresponds to i?6o (60 with respect to critical), and 
is substantially larger than R^oo and /?200- Cluster A^oop's are 
computed for th e S82 data using the max BCG prescription as 
implemented by Menanteau et al. ( 2010b| see Section 2.2 for 
details). Table |6jsummarizes the measured A^oop and inferred 
values of M500 and R^oo for each cluster field. 
A recent Planck study explores the SZE-optical scaling 



TABLE 6 
Optically Informed Properties 



Field 


N 200p 


M500 
(X10 14 M G ) 


-#500 

(Mpc) 


Y500 
(xl0~ 4 Mpc 2 ) 


ACTJ0022 
ACTJ2051 
A2631 


113± 11 
25 ±5 
77 ±9 


7.77 ± 1.12 
1.54 ±0.33 
5.19 ±0.75 


1.03 ±0.05 
0.73 ± 0.05 
1.11 ±0.05 


0.84I°;J! 
1.431^ 



Note. 



A^oop's are measured from SDSS S82 data, M SI 



is derived from A^oop using the scaling relation of Ro zo et aT] 
(20091, and Rsoo is inferred from the mass. This K500 is then 
used to compu te Ysoo from the SZA data for comparison to the 
Planck results (Planck Collaboration 201 lc I. 

relations by employing a multi-frequency matched filter on 
Planck maps at the positions of the SDSS maxBCG clusters 
( |Planck Col laboratio n |2011c| l. This work finds an offset be- 
tween the measured integrated Compton-y-optical richness 
relation compared to model predictions for the full maxBCG 
sample. However, when using the X-ray subsample, the au- 
thors find good a greement between the predicti on and the 



model. Following Planck Collaboration! (201 lc 1, we deter 
mine R 5Q o from the optical properties via the A^oop-^soo rela- 
tion (Table |6|l and compute F500 within that radius, using the 
fits of the ArO profile to the SZA data. Ysoo is rescaled to 
redshift z — using the evolution of the Hubble parameter 

for a flat universe, E(z) - ^Q.m(1 + z) 3 + Qa, and a fiducial 
distance D A = 500 Mpc, as 



^500 = Y 5 qoE{z) 



-2/3 



(Z) A (z)/500 Mpc) 2 [arcmin 2 ] 



(4) 



Figure |5]compares our two ACT clusters (red points) to those 
of Planck (black points) for the Fson-A^oop relation. The best 
fit power law for the Planck data over the full maxBCG cat- 
alog is also shown (line). The known X-ray cluster, A2631, 
lies above the Planck relation, though with large uncertainty. 

5. DISCUSSION 

The density of compact radio sources is higher in cluster re- 
gions compared to the fiel d (Coora y" gt al.|19 98 ; Mass ardi &| 
|De Zotti|2004HLin & Mohr|2007||Coble et al.|2007| l. Contam- 
ination of the SZE signal from radio sources could potentially 
bias flux and mass estimates from the SZE. The radio source 
fluxes extrapolated to 148 GHz (Table [4]i provide information 
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on potential contamination of the SZE decrement signal. The 
two sources with projected radius from the cluster R > 6' will 
have little impact on cluster detection and flux. ACT 148 GHz 
equatorial maps have a noise level of around 2 mJy beam -1 so 
that most sources are expected to fall well below the ACT 
map noise level, and the brightest source is expected to be 1 .5 
times the noise level. Very Large Array (VLA) observations 
of galaxies in nearby clusters between 5 and 40 GHz find that 
about 60% of the point source s show a fl attening of the spec- 
tral shape above 8 GHz ( |Lin et aL"1|2009| ). This implies that 
extrapolating from low frequency to high frequency yields a 
lower limit on the contaminating flux. A conservative upper 
limit on the radio source flux density at 148 GHz is obtained 
by increasing the extrapolated estimates by a factor of two. 
However, recent simulations of the microwave sky suggest 
that only 3% of clusters have their 148 GHz SZE dec rements 
contaminated at the > 20% level ( |Sehgal et al.|2010[ >. There 
is no indication in the ACT data of contamination by sources 
in these three cluster fields. The estimated 148 GHz flux den- 
sities suggest that radio sources do not significantly impact 
cluster detection in surveys for the brightest clusters, although 
they could potentially bias flux measurements in some frac- 
tion of the clusters at a level of 3-6 mJy (< 20% for typical 
-30 mJy integrated SZE fluxes). We reiterate that extrapo- 
lation to higher frequencies is uncertain and note that this is 
based on sources in only 3 cluster fields. More precise flux 
density estimates at ACT frequencies will be obtained through 
observing a larger number of ACT cluster fields over a range 

of frequencies. 

The VLA FIRST survey ([ Becker et aT1[T995] > covers the 
ACTJ0022 and A2631 fields! Three out of the 4 radio 
point sources have NVSS and FIRST flux densities that agree 
within 68% confidence. The 107 mJy NVSS source in 
A2631's field has a 91 mJy flux density in the FIRST cat- 
alog, resulting in a predicted 148 GHz flux density that is 
9% higher than that estimated from the NVSS flux density. 
These surveys are in different configurations and observe the 
same fields at different times, thus providing a rough handle 
on potential time-variability of these sources, which could im- 
pact the contamination of the SZE signal a nd ou r use of point 
sources to calibrate our SZA data (Section [3. 2| i. A 10% flux 
density variability is encompassed by our choice of a 20% 
calibration uncertainty for the SZA data. 

The S82 data for ACTJ2051 exhibit a red sequence typi- 
cal of galaxy clusters and provide a spectroscopic redshift. 
The optical mass estimates from the S82 data suggest that 
ACTJ2051 is less than a third of the mass of the other two 
clusters (Table |6j). Using the Planck T500--V200P relation, we 
estimate Ysoo = 2.7 + 1.2 x 10~ 6 Mpc 2 , which corresponds to 
415 ± 190 ,uJy of integrated SZE flux within R 500 at 31 GHz. 
This corresponds to an SZA signal smaller than 2<x for the ob- 
servations considered here, below the SZA detection thresh- 
old. Furthermore, the above significance calculation assumes 
that the entire SZE flux within ^500 is actually contained 
within the synthesized beam (effective resolution) of the SZA. 
However, the SZA synthesized beam is 1'.7, smaller than /?5oo 
(2'.6) for this cluster, and the cluster SZE signal will be diluted. 

Our mass estimates for ACTJ0022 and A2631 confirm that 
they are massive systems (M500 — 10 15 M Q ). This is con- 
sistent with optical and X-ray follow-up of the initial cluster 
results from ACT's southern region that sugg ests the high sig 



nificance detections will be > 8 x 10 14 M Q (Menanteau et al. 



20Tfy . The p and A10 model fits to the SZA data yield con- 
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Fig. 5. — Integrated Compton->'-optical richness relation, Fjoo-A^oOp. for 
Planck (black) and this work (red). Also shown is the best fit relation to the 
Planck data (solid line). Integrated Compton-j> values are spherical values 
within R500 scaled to redshift zero and a distance D k = 500 Mpc. Richness, 
A^oflp. i s the number of red sequence galaxies within R200p- 

sistent masses for both clusters. A comparison of the SZE, 
X-ray, and joint mass estimates of A263 1 shows more scat- 
ter. This scatter might be indicative of A263 1 being an unre- 
laxed cluster, especially in light of the asymmetric structure 
seen in both the SZA and Chandra observations (see Fig. [2]). 
However, given the uncertainties in the masses, the different 
methodologies yield consistent results for A263 1 as well. 

ACTJ0022 was discovered by ACT and we present ini- 
tial mass estimates but there are previous X-ray analyses of 
A2631. Temperatures are computed from regions that dif- 
fer between the different works. Using t he ear lier of the 
two Chandra observations, Maugh an et al.| (2008) finds T e — 
6.5 ± 0.6 keV within (0.15-l)fl 50 o, where R 50 o = L2 Mpc. 
Analyses of XMM-Newton data using a number of methods 
for computing Rsqq yields T e = 7.5 ± 0.5 keV (within (0.2- 
0.5)/?5oo) and an a verage M500 = (8 + 2) x 10 14 M Q within 
^500 = 1.2 Mpc (Zhang ~et al.||2008] ). Another analysis of 
XMM-Newton data finds T e = 7.7+0.6 keV within an aperture 
cont aining 90% of the ba ckground-subtracted surface inten- 
sity ( jAndersson et al.|2009| >. Our derived T e = 7.3 ± 0.6 keV 
agrees to within 1 .3cr with other Chandra and XMM- Newton 
measuremen ts, despite the differences in radii. The Zhang 
|et al.| ( |2008| l M500 is consistent with our derived SZE and X- 
ray masses to within 2<x in the most discrepant case. Despite 
data from various observatories and different methodologies, 
the temperatures and masses of A263 1 are consistent. 

A number of systematics may affect these mass measure- 
ments, summarized in Table [7] as percentages for a single 
cluster. Both simulation and analysis of X-ray data suggest 
that asphericity typically affects mass estimates at the 5-10% 
level when measured at large radii such as ^500 (e.g., Math- 
iesen et aL]|1999| [PTffaretti et al. 2003). We conservatively 
adopt 10% for the effects of asphericity. Simulations sug- 
gest that small-scale fluctuations, often called clumping, wil l 
cause a 10% overestimate of mass (Mathiesen et al. 1999). 
Non-thermal pressure support may affect both SZE and X- 
ray mass estimates, causing mass to be underestimated at the 
10-20% level (e.g., |Nagai et al.||2007| [Lau et al.||2009l >. For 
the SZE-only analysis, changes in the assumed gas fraction, 
on average, change the mass by 12%, with lower gas fraction 
leading to high er mass; the fixed metallicity assumption is a 
< 1% effect ( |Mroczk owski 2011), which we neglect. Rela- 
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TABLE 7 
Potential Systematics on Mass 



Source 


SZE 

(%) 


X-ray 

(%) 


Calibration 


±10 


±10 


Asphericity 


±10 


±10 


Clumping 


-10 


-10 


Non-thermal P e 


+ 15 


+ 15 


Assumed f g 


±12 




Relativisic corrections 


+3 




Total" 


+24 
-21 


+21 
-17 



■ Added in quadrature. 

tivistic corrections to the SZE are a 3% correcti on at 30 GHz 
for an 8 keV cluster like A2631 ( |Itoh et al.|1998) . 

The SZA and Chandra calibrations both affect the mass es- 
timates. The SZE mass determinations depend on the SZE 
signal asMcc (AT) 1 / 2 so that the conservative 20% calibra- 
tion results in a 10% change in mass. Recent changes in the 
Chandra calibration can change the cluster temperatures in- 



ray 



ferred from spectroscop ic fits by 10% ( Reese et al.| | 2010 
|Nevalainen et al.|2010 i. This will directly impact our X 
mass estimates since M oc T e , 

Potential systematics are summarized in Table[7] with totals 
added in quadrature. Both SZE and X-ray estimates have po- 
tential systematics at the 20% level, roughly the same order 
as the statistical uncertainties (Table |5j. 

6. CONCLUSION 

We obtained SZA follow-up observations of three optically- 
confirmed clusters from preliminary ACT maps along the 
celestial equator. ACT-CL J0022-0036 is a massive, high- 
redshift cluster newly discovered by ACT. The SZA detects 
two of the three clusters at high significance. ACT and SZA 
data show good qualitative agreement (Figure fT). The clus- 
ter A2631 shows good agreement between SZE and X-ray 
data (Figure [2]) with both peaks aligning and similar E-W 
elongation seen in both wavebands. Initial mass estimates 
confirm that ACTJ0022 and A263 1 are massive clusters with 
M 5 oo - 10 15 M Q . 

Two compact radio sources are detected by the SZA in each 
cluster field (Table 0}. Using NVSS 1.4 GHz flux densities, 
we compute spectralindices and predict the flux densities in 
ACT's 148 GHz band. The radio sources are expected to be 
< 6 mJy in the ACT 148 GHz maps, suggesting that radio 
sources are not a significant contaminant for detection of high 
mass clusters. However, they can still impact the measured 
SZE signal. As an example, the brighter source in A2631 
could be filling in the SZE signal at the < 20% level, assuming 
the extrapolation of the sources' lower-frequency flux densi- 
ties holds. A more precise determination of potential contam- 
ination of the SZE signal from compact sources at ACT fre- 
quencies will be obtained through observing a larger number 
of ACT cluster fields over a range of frequencies. 

Optimized for different purposes, ACT and SZA provide 
complementary SZE data on galaxy clusters. With the ability 
to quickly image large regions of the sky, ACT is well-suited 
to finding clusters. Targeted observations with the SZA al- 
low deep integrations for detailed cluster studies. The SZA 
images (Fig. [j} show peak S/N of 12 and 10 for ACTJ0022 
and A2631, respectively, compared to S/N 6 for the match- 
filtered ACT images. The SZA provides higher resolution 
than ACT, measuring angular scales 15"— 5', well-matched to 



clusters. In addition, the spatial filtering of the interferome- 
ter provides a clean method of disentangling cluster emission 
from radio point source emission. 

This pilot study of SZA follow-up observations of ACT- 
detected clusters shows that the detected clusters are massive 
systems. Cluster abundances are e xponentially sensitive to 
mass (e.g., |Press & Schechter||1974| l, with the most massive 
clusters providing the most leverage on cosmological studies. 
Finding extremely massive systems at high redshif t can poten- 
tially rule o ut the current A-CDM paradigm (e.g., Mortonson 
et al.|20lT I. In addition, cosmological determinations utiliz- 
mg massive systems are expected to be less susceptible to the 
effects of non-gravitational astrophysics. The highest mass 
systems from ACT's two survey regions will comprise our 
core sample for derivation of cosmological parameters from 
clusters. It is crucial to measure well the integrated SZE sig- 
nal and the masses of these clusters for proper cosmological 
interpretation of cluster yields from current and future SZE 
cluster surveys. Deep, targeted SZA observations provide a 
method of disentangling the point source and cluster emis- 
sions and enable robust estimates of the integrated SZE signal, 
7500, and initial mass estimates. When combined with mass 
measures from X-ray and weak lensing, this will provide a 
robust measure of the mass-SZE flux scaling relation. 
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